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1. Introduction 

As it is well known, general relativity theory (GR) does not lead to restrictions for 
admissible energy densities in the case of gravitating systems with positive values of 
energy density and pressure satisfying energy dominance condition. As a result, the 
problem of gravitational singularities takes place in the frame of GR [1]. The appearance 
of nonphysical state with divergent energy density limiting world lines in the past or 
in the future is characteristic feature of Friedmannian cosmological models (problem 
of cosmological singularity - PCS) and also of collapsing systems. From physical 
point of view, the problem of gravitational singularities is connected with the fact 
that gravitational interaction in the case of gravitating systems with positive energy 
densities and pressures in the frame of GR as well as Newton’s theory of gravitation 
has the character of attraction but not repulsion. Note that gravitational interaction in 
GR can have the repulsion character in the case of gravitating systems with negative 
pressure (for example systems including massive or nonlinear scalar fields). According 
to accepted opinion, it is possible cause of acceleration of cosmological expansion at 
present epoch. However, such effect does not permit to solve the PCS in the frame of 
GR [2]: all Friedmannian cosmological models of flat and open type, and the most part 
of closed models are singular. 

There were many attempts to solve the problem of gravitational singularities and 
at first of all the PCS in the frame of GR and other classical theories of gravitation; a 
number of particular regular cosmological solutions was obtained (see [3] and references 
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given herein). In connection with this note that the solution of PCS means not only 
obtaining regular cosmological solutions, but also excluding singular solutions; as a result 
regular behaviour of cosmological solutions has to be their characteristic feature. The 
most existent attempts to solve the PCS do not satisfy indicated conditions. According 
to wide known opinion, the solution of PCS and generally of the problem of gravitational 
singularities of GR has to be connected with quantum gravitational effects, which must 
be essential at Planckian conditions, when energy density is comparable with Planckian 
one. Previously some regular bouncing cosmological solutions were obtained in the 
frame of candidates to quantum gravitation theory — string theory/M-theory and loop 
quantum gravity (see, for example, [4-6]). From physical point of view, these solutions 
have some difficulties [7]. In the case of bouncing cosmological solutions built in the 
frame of string theory the condition of energy density positivity for gravitating matter 
is violated. In the case of loop quantum cosmology a bounce takes place for microscopic 
model having a volume comparable with Planckian one. If one supposes that Universe 
at compression stage is macroscopic object, one has to explain the transformation of 
macro-universe into micro-universe before a bounce, this means one has to introduce 
some model inverse to inflation. 

As it was shown in a number of papers (see [7, 16] and references given herein), 
gauge theories of gravitation (GTG) and at first of all the Poincare GTG, which are 
natural generalization of GR by applying the gauge approach to gravitational interaction 
[8,9], offer an opportunity to solve the PCS. All homogeneous isotropic cosmological 
models including inflationary models are regular in metrics, Hubble parameter, its time 
derivative by certain restriction on equation of state for gravitating matter at extreme 
conditions. Unlike Friedmannian models of GR nonphysical state with divergent energy 
density does not appear because of gravitational repulsion effect at extreme conditions, 
which takes place in GTG in the case of usual gravitating systems with positive energy 
densities satisfying energy dominance condition. 

The present paper is devoted to analysis of gravitational repulsion effect at extreme 
conditions in the frame of GTG in the case of homogeneous isotropic gravitating systems. 
In Section 2 some important properties of homogeneous isotropic models in the frame of 
Poincare GTG are briefly discussed. In Section 3 generalized cosmological Friedmann 
equations without specific solutions for such models filled by usual matter and scalar 
fields are deduced. In Section 4 extreme conditions leading to gravitational repulsion 
effect are analyzed. 

2. Homogeneous isotropic gravitating models in Poincare GTG 

In the framework of gauge approach to gravitation the Poincare GTG is of the greatest 
interest [8,7]. Gravitational field variables in Poincare GTG are the tetrad (translational 
gauge field) and Lorentz connection (Lorentz gauge field), corresponding field strengths 
are torsion and curvature tensors, and physical space-time is 4-dimensional Riemann- 
Cartan continuum. As sources of gravitational field in Poincare GTG are covariant 
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generalizations of energy-momentum and spin tensors. Unlike gauge Yang-Mills fields, 
for which the Lagrangian is quadratic in the gauge field strengths, gravitational 
Lagrangian of Poincare GTG can include also linear in curvature term (scalar curvature), 
which is necessary to satisfy the correspondence principle with GR. By using sufficiently 
general expression for gravitational Lagrangian, homogeneous isotropic gravitating 
models were investigated in the frame of Poincare GTG in a number papers (see for 
example [7, 10, 16]). Because of high spatial symmetry gravitational equations depend 
weakly on the structure of quadratic part of gravitational Lagrangian that permits to 
obtain conclusions of general character. Below some important relations describing 
homogeneous isotropic models in Poincare GTG will be given. 

Spatial symmetries of Riemann-Cartan space are defined by a set of linearly 
independent Killing’s vectors, with respect to which the Lie derivatives of metric and 
torsion tensors vanish (see [1] Chapter 2, [17]). Homogeneous isotropic models possess 6 
linearly independent Killing’s vectors and metric tensor in co-moving system of reference 
has the form of Robertson-Walker metrics [17]: 

9nv = diag ^ 1, -, - R 2 (t) r 2 , - R 2 (t) r 2 sin 2 0 ^ , 


where R(t) is the scale factor, k = +1,0,—1 for closed, flat, open models respectively 
(the light velocity c = 1), spatial spherical coordinates are used. Then the torsion tensor 
S x pi/ = —S x up satisfying symmetry conditions can have the following non-vanishing 
components [10]: S 1 10 = S 2 20 = S 3 30 = S(t ), S 123 = S 231 = S 3i2 = S(t ) ^/==f sing, 
where S(t) and S(t) are functions of time. The functions S and S have different 
properties with respect to transformations of spatial inversions, namely, the function 
S(t) has pseudoscalar character. By supposing S = 0, we obtain for the curvature 
tensor F^ v ap = —F v ^ ap = —F^ u pa the following non-vanishing components: K 01 oi = 
F 02 o 2 = F 03 os = A and F 12 i 2 = F 13 13 = F 23 23 = B with 


A = 


R-2RS 


B = 


k + (R - 2 RS 


R R 2 ’ 

and Bianchi identities in this case are reduced to the only relation 


( 1 ) 


B + 2H(B -A) + 4AS = 0, (2) 

where H — is the Hubble parameter, and a dot denotes differentiation with respect 
to time. 

We will use gravitational Lagrangian in the form 


C G = f 0 F + F a ^ (/, F aPpv + f 2 F apPu + f 3 F pua p) + (/ 4 F pv + / 5 F vp ) + / 6 F 2 

+S a ^ ( ai S apu + a 2 S vpa ) + a 3 FVS/" 3 , 


where F pv = F a pa „, F = f\ (i = 1,2, ...,6), a k (k = 1,2,3) are indefinite 
parameters, /o = (167rG) -1 , G is Newton’s gravitational constant. Gravitational 
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equations for homogeneous isotropic gravitating models (with S — 0) are reduced to 
the system of 3 equations [10], which by using (2) can be written as 


6 f 0 B - 12/ (A 2 - B 2 ) -3 a(H-S)S = p, 

2/o (2A + B) + 4/ (A 2 - B 2 )-a(s + HS - S 2 ) = -p, ( 3 ) 

f (A + B^j + [/ 0 + |a + 4/ (A + B)~\ S = 0. 

where / = ,f\ + |/ 2 + h + .A + h + 3/ 6 , a = 2a 1 + a 2 + 3a 3 , p is the energy density, p 
is the pressure and the average of spin distribution of gravitating matter is supposed to 
be equal to zero. The system of equations (3) leads to cosmological equations without 
high derivatives if a = 0 [10] (see below). Then we find from (3) the curvature functions 
A and B and torsion S in the following form 


A = _ P + 3p - ft (p - 3 pf /2 

12/ 0 1 + a (p - 3p) 

B = 1 p + a(p-3p) 2 /4 

6/o 1 + a (p - 3 p) 

S(t) = -~ In \l + a(p - 3p)\ , 


where indefinite parameter a = 


/ 


3/o : 


has inverse dimension of energy density. Note, that 


in the case S A 0 Bianchi identities for homogeneous isotropic models are reduced to two 
relations and system of gravitational equations of Poincare GTG includes 4 equations. 
However, these equations have solution (4) together with S — 0 always, if a = 0 [18]. 


3. Generalized cosmological Friedmann equations in GTG 


By using expressions (1) of curvature functions for homogeneous isotropic gravitating 
models and the solution (4) of gravitational equations of Poincare GTG we obtain the 
following generalized cosmological Friedmann equations (GCFE) 


R- 1 ^- 


d 

dt 


k 

R 2 


dR „ d 
— + R— 

dt dt 


d_ 

dt 


In 


R\f 11 + a (p - 3p)| 


8? tG P + f (p ~ 3p) 2 
3 1 + a (p — 3 p) 


In V 11 + a (p — 3p)j 


47tG p + 3p - f (p - 3p) 


1 + a (p — 3p) 


(5) 


( 6 ) 


The difference of (5)-(6) from Friedmannian cosmological equations of GR is connected 
with terms containing the parameter a. These terms arise from quadratic in 
the curvature tensor part of gravitational Lagrangian, which unlike metric theories 
of gravitation does not lead to high derivatives in cosmological equations.! The 
conservation law in GTG has usual form 


p + 3H (p + p) = 0. (7) 

| As it was shown in [11,12] generalized cosmological Friedmann equations (5)-(6)can be deduced also 
in the frame of the most general GTG - affine-metric GTG [9]. 
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Now by using GCFE (5)-(6) we will consider homogeneous isotropic models filled 
by non-interacting scalar field f minimally coupled with gravitation and gravitating 
matter with equation of state in general form p m = p m (p m )- (The generalization for the 
case with several scalar fields can be made directly). Then the energy density p and 
pressure p take the form 


P = ^4> 2 + V + p m (p > 0), p = ^f 2 - V + p m , 


( 8 ) 


where V = V {f) is a scalar field potential. By using the scalar field equation in 
homogeneous isotropic space 


dV 

* +3H * = ~ai 

we obtain from (7)-(9) the conservation law for gravitating matter 
Pm T 3 H ( p m Pm) = 0. 

By using (8)-(10) the GCFE (5)-(6) can be transformed to the following form 


(9) 

( 10 ) 


H 


Z - On ( e/ + -Y 


3V ■ 

+ 3 “^ 


. k y- 

+ w- 


8ttG 


Pm + ~f 2 + V + -a ( 4V — (j) 2 + p r 


3 Pr 


Z, 


( 11 ) 


H 


Z + 3a \ 4> 2 + -Y 
—3aY 


-12 aH 


Z + H' 2 yZ - 15 af) 2 
9 a f dp , 


-18a 1 


2 

dV • 

dV 

df 


- - Y + 3 {p m + Pm) 2 ~Y~ 2 

dpm dp r 


Z - 18a 2 


+ 2 Y 


Z + 3a ( + —Y 


3 a 


d 2 V ■„ 


df> 2 


dV_ 

df 


8i tG 


Y f ^ {Pm T 3 Pm) 


1 / . \2 
+ -a {4V - f 2 + p m - 3p. 


( 12 ) 


+ pm- 3 Pm ) and Y = (p m + p m ) ( 3 - 1 


where Z = 1 + a \ 4V — a2 

By transformation of GCFE (5)-(6) to the form (11)-(12) these equations were 
multiplied by Z. As a result (11)-(12) have specific solutions when Z — 0 [13, 7]. In 
fact, by using the expression of the Hubble parameter following from (11) 


H + = 


dV ■ r~ 
-3a —— f ± vD 
df 


( 13 ) 


where 


D = 


8ttG 


Z + 3a[f z + -Y 


Pm + ~& 2 + Y + —a ( 4V — f 2 + p n 


3 p r 


7 ^72 

z ~lv z ’ 
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it is easy to show that (12) is satisfied, if Z = 0. ft is because the terms in (12), which 
don’t include Z as multiplier, vanish by virtue of (13). Excluding these terms and 
dividing obtained equation on Z we will have, instead of (12), the following equation 
without specific solutions 


H 


Z + 3a ( <f>~ + — Y 

3a / dp. 


+ 3 H 2 


2 

= 87 tG 


1 -Y + 3 ( p m + p m ) 
dpm d,p 2 m 

1 , ,1 


+ olY 
2 d 2 p m 


3a 


d 2 V 


d<f 2 


dV 

dcj) 


V + — (p m — Pm) + ~ O' — (j) 2 + p rn — 3/J m J 


2k 

B? Z ' 


(14) 


The system of equations (11) and (14) is equivalent to GCFE (5)-(6) in considered case 
of models filled by usual matter and scalar fields. If the interaction between scalar 
holds and usual gravitating matter is not neglecting, is allowed the GCFE (11), (14) 
have to be generalized. By taking into account the interaction by means of scalar held 
potentials, which in this case depend also on the energy density of gravitating matter 
p rn , namely V = V (0, p rn ) [7], we can obtain the generalization of equations (11), (14) 
by similar way. The form of (11) does not change, but the value of Y in this case is 
defined as Y = (3p m — 417) — lj. The equation (11) is generalized as 


H 


Z + 3a U 2 + iy) 

+ 3 H 2 

V 2 ) 



+ —Y 
2 


3 a Y 


1 + 


dPm , Pm + P m d 2 V\ 3a (p m + p m ) 2 d 2 


+ 


2 1 + £ V dpm l + d^: d P 2 m 


d 2 V 




(3 p m - 417) 


27 Tj 1 Pm "h Pm 

D — 2 oT O 

2 / . dvY d(j)dp n 


1 + 


8pn 


1 d , 


l>17 

fdV\ 2 ' 

1 + o , (Pm + 217) 

S Cl Prn 

~\~ 3a 

1 

Cv, ' 
-S- . 

to 

"O- 

(n 

1 

( 90 ) _ 


= 87 tG 


1 1 / N 2 

V + - {pm - Pm) \ - 4> 2 + Pm. ~ 3 p, 


2k 

R 2Z ' 


(15) 


The system of (11), (15) does not have specihc solutions (with Z = 0) unlike (9)-(10) 
of reference [7]. By using obtained equations (11), (14), (15) the repulsion gravitational 
effect can be analyzed in the case of homogeneous isotropic gravitating systems in the 
frame of GTG. 


4. Gravitational repulsion effect in GTG 

As it was noted above, the difference of GCFE from Friedmannian cosmological 
equations of GR is connected with terms containing the parameter a. The value of 
|a| _1 determines the scale of extremely high energy densities. Solutions of GCFE (5)- 
(6) coincide practically with corresponding solutions of GR if the energy density is 
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small | a(p — 3p)\ <C 1 (p f \p). The difference between GR and GTG can be essential 
at extremely high energy densities \a(p — 3p)\ > 1. Ultrarclativistic matter ( p = |p) 
and gravitating vacuum (p = — p) with constant energy density are two exceptional 
systems because GCFE (5)-(6) are identical to Friedmannian cosmological equations of 
GR in these cases independently on values of energy density, and non-einsteinian space- 
time characteristics vanish. Properties of solutions of equations (5)-(6) at extreme 
conditions depend on the sign of parameter a and certain restriction on equation of 
state of gravitating matter. The study of inflationary models including scalar fields 
shows that GCFE (5)-(6)lead to acceptable restriction for scalar field variables if a > 0 
[7, 16]. In the case a > 0 all cosmological solutions including inflationary solutions have 
regular bouncing character if at extreme conditions p > |p. There are physical reasons to 
assume that the restriction p > is valid for gravitating matter at extreme conditions 
[14], Really in the case of perfect gas of fermions at zeroth absolute temperature, the 
pressure tends to the value p/3 from below, if density of gas increases. Then we have 
for nuclear matter at extreme conditions p > \p because of strong nucleon interaction 
[14], We will suppose below that the condition p > is valid for gravitating matter 
at extreme conditions and, in particular, at the beginning of cosmological expansion. § 
Note, that this condition is valid for so-called stiff equation of state p = p used in the 
theory of the early Universe (Ya. B. Zeldovich and others). 

The GCFE lead to restrictions on admissible values of energy density. In fact, if 
energy density p is positive and a > 0, from equation (5) in the case k = +1, 0 follows 
the relation: 

Z = 1 + a (p — 3 p) > 0. (16) 

The condition (16) is valid not only for closed and flat models, but also for cosmological 
models of open type (k = —1) [7]. In the case of models filled by usual gravitating 
matter without scalar fields the equation Z = 0 determines limiting (maximum) energy 
density, and regular transition from compression to expansion (bounce) takes place for all 
cosmological solutions by reaching limiting energy density. Near limiting energy density 
the gravitational interaction has the repulsion character. In the case of systems including 
also scalar fields a bounce takes place in points of so-called ’’bounce surfaces” in space 
of variables (0, f>, p m ) [7]. Near bounce surfaces as well as bounds Z = 0 gravitational 
interaction has the character of repulsion, but not attraction. By using GCFE in the 
form (11), (14) we will find below, by what conditions gravitational repulsion effect 
takes place. 

a) At first we will consider gravitating systems filled by gravitating matter without 
scalar fields. Then the acceleration a = = H + H 2 from (13)-(14) takes the following 

§ In the case a < 0 GCFE (5) - (6) lead to regular bouncing solutions if at extreme conditions p < 

[10, 3]. However, in this case GCFE have also singular solutions for some hypothetical superdense 
systems [15, 3]. 
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form 


a — ( z + -ay 


-3 


4v tG 


Pm Pm T 2 (Pm 3p r 


-2D 


7 , 3 v 9 « 
Z + -ar-— 

4 4 


dpn 


+ 3(j> m + p m ) 


_ 2/c _ 
R 2 Z 

2 d 2 p, 


Z + -aY 


^ i r m/ to 

Upm dpm, 

Obviously the repulsion (a > 0) will take place, if the expression in figured parentheses in 
(17) is positive. The domain of admissible energy densities, by which the repulsion effect 
takes place, depends on equation of state of gravitating matter at extreme conditions 
and on the value of parameter a. In particular, in the case of flat models (k = 0) with 
linear equation of state p rn = wp m (w = const > |) we obtain from (17) the following 
condition for energy densities corresponding to repulsion effect 


(17) 


y(x) = -(9 w 2 — l)x 3 + 


3(3tc 2 — 1) 2 3(9w + 5) 


x + 


-x - 1 > 0, 


(18) 


3w + 1 ' 2(3 w + 1) 

where x = ap m (3w — l) > 0. Moreover from condition Z = l—x > 0 follows, that x < 1. 
By this restriction the cubic equation y(x) = 0 has the only real root Xq and inequality 
(18) is fulfilled at x > x 0 . Numerical solution of equation y(x) = 0 gives the dependence 
of value xq on parameter w (see Fig.l). The gravitational repulsion effect in considered 
case takes place at energy densities defined by the following condition x 0 < x < 1, and 
the value x — 1 corresponds to limiting energy density. 


Xo 



Figure 1. The root xq as function of parameter w. 


b) In the case of systems including non-interacting gravitating matter and scalar 
held the condition of gravitational repulsion obtained from (11) and (14) is 
1 


87rG 


^ + 2 (P r 


a 


Pm) + j{ 41/ 


0 H“ Pm 3 Pm) 
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-3a 


dV ■ 


fdV 


- 


\d(j) 



Zi + 3of I + — Y 


2 r 


-3a— 


, 9« f dp m , x 

+ I ^ + 3 (Pm + Pm) , 2 
2 \ dp m dp n 


—2Z + 6a0 2 — ^aY' 
2 d 2 p n 


2 
> 0 


(19) 


Inequality (19) together with condition Z > 0 determine the domain of variables (j) 
and p rn at extreme conditions near bounds (Z = 0) and bounce surfaces in space of these 
variables, where gravitational repulsion effect appears. This domain is different for H + - 
and //--solutions of (11) and (14) corresponding to two values of the Hubble parameter 
(13) [7]. Note that (19) describes also repulsion effect at small energy densities like to 
GR, when the pressure is negative because of contribution of scalar held. In the case 
of systems filled by interacting gravitating matter and scalar held the generalization of 
condition (19) can be obtained from (11) and (15). Note, that in the case of closed and 
open models conditions for gravitational repulsion effect include the term with scale 
factor R depending on 3-space curvature also, this means that the repulsion effect for 
homogeneous isotropic systems depends on global structure of gravitating model. 


5. Conclusion 

The analysis of gravitational repulsion effect at extreme conditions in GTG presented 
above shows that this effect takes place near limiting energy density, near a bounce, 
and it depends on content and properties of gravitating matter at extreme conditions 
(equation of state for gravitating matter, the form of scalar held potentials etc) and also 
on the scale of extremely high energy densities defined by parameter a. If the value of 
limiting energy density is essentially less than the Planckian one, gravitational repulsion 
effect appears at classical conditions, when the quantum gravitational corrections 
according to widely known opinion are not essential. || In the case of systems including 
scalar helds the value of limiting energy density is different for different solutions and can 
be essentially greater than a -1 , but the appearance of gravitational repulsion effect does 
not depend on this fact [16]. If limiting energy density is comparable with the Planckian 
one, quantum gravitational corrections have to be taken into account, although classical 
GTG ensures satisfactory non-singular behaviour of gravitating systems. Because GCFE 
in the case of gravitating systems with small energy densities lead to consequences 
similar to that of GR, gravitational repulsion effect at such conditions can appear in the 
case of gravitating models with negative pressure (dark energy) like to GR. In particular, 
if GCFE (5)-(6) include cosmological constant, the value of which corresponds to dark 
energy density at present epoch, regular cosmological solutions of (5)-(6) will contain 
the stage of accelerating cosmological expansion by dominating of vacuum energy 

|| Probably, in the frame of considered GTG the quantum properties of gravity can be important near 
the bound Z — 0, where the values of the torsion function S are extremely large. 
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(cosmological constant). 

The study of gravitating systems in the frame of GTG shows, that important 
features of gravitational interaction at extreme conditions depend essentially on 
properties of gravitating matter, which are determined by other fundamental physical 
interactions. To describe the evolution of gravitating models we have to know how 
corresponding properties of gravitating matter (at first of all its equation of state) change 
by model evolution. This conclusion obtained by investigation of homogeneous isotropic 
models has sufficiently general character, because the form of used GCFE (5)-(6) does 
not depend on detailed structure of quadratic in the curvature part of gravitational 
Lagrangian of GTG. At the same time the dynamics of gravitating systems depends on 
the structure of quadratic part of gravitational Lagrangian of GTG, if the homogeneity 
and isotropy are broken. In connection with this, the search of GTG leading to the 
most satisfactory consequences in general case of inhomogeneous anisotropic models is 
of direct physical interest. 
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